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Higher order Laguerre-Gauss (LG) beams have been proposed for use in future gravitational wave 
detectors, such as upgrades to the Advanced LIGO detectors and the Einstein Telescope, for their 
potential to reduce the effects of the thermal noise of the test masses. This paper details the theo- 
retical analysis and simulation work carried out to investigate the behaviour of LG beams in realistic 
optical setups, in particular the coupling between different LG modes in a linear cavity. We present 
a new analytical approximation to compute the coupling between modes, using Zernike polynomials 
to describe mirror surface distortions. We apply this method in a study of the behaviour of the 
LG33 mode within realistic arm cavities, using measured mirror surface maps from the Advanced 
LIGO project. We show mode distortions that can be expected to arise due to the degeneracy of 
higher order spatial modes within such cavities and relate this to the theoretical analysis. Finally 
we identify the mirror distortions which cause significant coupling from the LG33 mode into other 
order 9 modes and derive requirements for the mirror surfaces. 

PACS numbers: 04.80.Nn, 95.75.Kk, 07.60.Ly, 42.25.Bs 
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I. INTRODUCTION 

The sensitivities of second generation gravitational 
wave detectors such as Advanced LIGO and Advanced 
Virgo are expected to be limited by the thermal noise 
of the test masses within a significant range of signal fre- 
quencies around 100 Hz [1 . To reach even better sensitiv- 
ities, it has been proposed to use laser beams with an in- 
tensity pattern other than that of the fundamental Gaus- 
sian beam to reduce the effects of this thermal noise l2||3j. 
A beam whose intensity is distributed more homoge- 
neously over the mirror surface, for the same clipping 
losses, benefits from a more effective averaging over the 
mirror surface distortions caused by thermal effects [4]. 
The specific advantage of using higher order LG modes, 
as opposed to mesa [5 and conical [6] beams, is that they 
are compatible with spherical mirrors as currently used 
in GW detectors and other high precision optical setups. 
Research into the potential of the LG33 mode in gravita- 
tional wave detectors has been carried out using numer- 
ical simulations and table-top experiments [7^, ^. The 
sensing and control signals for an LG33 beam were found 
to perform as well as for the fundamental mode in all 
aspects examined and the LG33 behaved as expected in 
short linear and triangular optical cavities. However, an 
optical cavity resonant for a higher-order Gaussian mode 
is degenerate so that a number of modes can resonate at 
the same time. This is a fundamental difference to a well 
designed cavity for the fundamental Gaussian mode, in 
which any resonant enhancement of other modes can be 
suppressed. This degeneracy can potentially cause addi- 
tional optical losses. Simulations have shown that the use 
of the LG33 beam, compared to the fundamental mode, 
LGoo, could result in a significant contrast defect at the 
dark fringe [9h11j. It is the aim of this paper to inves- 
tigate how mirror surface distortions affect the purity of 
an LG33 beam in high finesse cavities, by analytical cal- 
culation and numerical simulation. We will focus on the 



direct coupling from a distorted mirror and how this af- 
fects the mode content in a linear cavity, with our final 
aim to produce specifications for the mirror surfaces. 



II. REPRESENTING MIRROR SURFACE 
DISTORTIONS AND LG MODES 

A. Mirror surface maps 

Ideally the mirrors in gravitational wave detectors 
should be perfectly smooth with a radius of curvature 
matching that of the incident beam. However, real mir- 
rors deviate from a perfect surface, altering the beams 
which interact with them. If a beam U{x^y^ z) is incident 
on a distorted surface described by Z(x, y) and uniform 
reflectivity r, then the reflected beam is given by: 

Uref{x,y,z) = U{x,y,z) rexp{i2kZ{x,y)), (1) 

Fig. [1] illustrates this effect. 

In order to investigate the effects of surface distortions 
measured mirror surface maps can be used. The term 
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Figure 1: A diagram illustrating the phase shift introduced 
when a beam is reflected from a distorted mirror surface with 
reflectivity r. The surface Z{x,y) is defined as a height field 
across a plane perpendicular to the optical axis along z. 



mirror map refers to an array of data detailing the op- 
tical properties of a mirror, often its surface height in 
nanometers. This data can be used to represent realistic 
mirrors in numerical simulations of gravitational wave de- 
tectors. Mirror maps have been produced from uncoated 
Advanced LIGO mirror substrates which represent the 
best mirror surfaces of this kind currently available. In 
the following we have made use of one such map, the sur- 
face map of the substrate ETM08 [12 . The deviation of 
this map surface from a perfectly spherical surface with 
radius of curvature 2249.28 m is shown in Fig. [2j This 
substrate shows an RMS surface figure error of 0.523 nm. 



B. Zernike polynomials 
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Zernike polynomials are well suited for the purposes 
of describing mirror surface distortions. Zernike polyno- 
mials can be used to describe classical distortions such 
as tilts and curvatures [13]. They are a complete set of 
functions which are orthogonal over the unit disc and de- 
fined by radial index, n, and azimuthal index, m, with 
n > m > 0. For any index m we have one odd and one 
even polynomial [14 : 



Z+'^(p, 0) = A^"^ cos{m(j))K^{p) even polynomial 
Z~'^(p, 0) = A""^ s\n{m<j))K^{p) odd polynomial 



(2) 



where p is the normalised radial coordinate, (j) is the az- 
imuthal angle, A^^ is the amplitude and R^{p) is the 
radial function. The radial function is given by the fol- 
lowing sum: 
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(3) 
for n — 771 even and otherwise. This gives n + 1 non- 
zero Zernike polynomials for each value of n (for tti = 
the odd polynomial is zero). Fig. [3] shows the surfaces 
described by the Zernike polynomials corresponding to 
orders (n) to 4. The lower order polynomials repre- 
sent some common optical distortions, some of which are 
summarised in table [H 

The odd polynomial describes the same surface as the 
even polynomial, but rotated by 90°. Combinations of 



Common name 



Offset 

1 zbl Tilt in x/y direction 

2 Curvature 

2 ±2 Astigmatism 

3 zbl Coma along x/y axis 



Table I: Summary of some common names for the lower order 
Zernike polynomials [13] . 



Figure 2: A surface plot of the mirror map ETM08 corre- 
sponding to the surface heights of an Advanced LIGO end 
test mass. The curvature has been fitted and removed from 
the map data. 



the odd and even polynomials relate to this same dis- 
tortion rotated by a given angle. The magnitude of this 
surface distortion can be given by the root mean squared 
amplitude of the polynomials: 



{A-„ 



{At"-- 



(4) 



where + refers to the even polynomial and — refers to 
the odd polynomial. Any surface defined over a disc can 
be described by a sum of Zernike polynomials, with the 
higher order polynomials representing the higher spatial 
frequencies present in the surface. 



C. Laguerre-Gauss modes 

The shape of any paraxial beam can be described as 
a sum of Hermite-Gauss or Laguerre-Gauss modes. The 
Laguerre-Gauss modes are a complete and orthogonal set 
of functions defined by radial index p and azimuthal in- 
dex /. The helical type of LG modes are typically given 
as [15 : 
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where k is the wavenumber, w{z) is the beam spot size 
parameter, ^{z) is the Gouy phase and Rc{z) is the ra- 
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Figure 3: Plots of the non-zero Zernike polynomials from n — 
to n = 4 with the odd polynomials with m = —n on the far 
left and the even polynomials with ?7i = n on the far right, 
in steps of 2. The colour scale represents negative surface 
heights with greens and blues, zero with black and positive 
surface heights with reds and purples. 



dius of curvature of the beam. Lp {x) refer to the as- 
sociated Laguerre polynomials. When considering these 
beams in cavities we note that the resonance conditions 
for these beams differ from that of a plane wave due to 
the (2p + |/| + l)'^{z) phase shift. The order of an LG 
mode is given by 2p+ |/| and modes with the same order 
will acquire the same round trip phase shift whilst circu- 
lating in a cavity. Therefore the cavity is degenerate for 
LG modes of the same order. 

The effects of mirror surface distortions on the shape 
of a reflected beam can be described in terms of coupling 
between LG modes. When a perfectly aligned Gaussian 
beam is reflected by a perfectly spherical mirror with the 
radius of curvature of the mirror matching that of the 
beam's phase front, the shape of the reflected beam is 
identical to the shape of the incident beam, or, in other 
words, the mode composition has not changed. However, 
if the mirror surface is distorted the reflected beam will 
generally have a different mode composition. The cou- 
pling from an incident mode (indices p and /) impinging 
on a completely reflecting surface Z into a mode (indices 
p' and V) in the reflected beam can be described by a 
coupling coefficient [161 Ull • 






f Up,iexp{2ikZ{r,^))U;.^, (6) 
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Z describes the surface height of the mirror and S de- 



scribes an infinite plane perpendicular to the optical axis. 
Currently the fundamental mode, LGqo, is used in 
gravitational wave detectors. Investigations have shown 
that mirror surface distortions have little effect on the 
beam purity when LGqo is used. The presence of mirror 
surface distortions introduces modes into the detectors 
other than the input mode, but since LGqo is the only 
mode of order any coupling out this mode will result in 
modes of a different order, which will be suppressed in the 
cavities. The LG33 mode is one of several modes of order 
9. In total there are 10 order 9 modes; LGo,±9, LGi,±7, 
LG2,±5, LG3,±3 and LG4,±i (Fig. |4|. These modes will 
be resonant in the arm cavities of the detectors, poten- 
tially resulting in a large proportion of the circulating 
power being in modes other than LG33. The distortions 
present in the mirrors in each arm cavity will be different 
and so the mode content in each arm will differ, resulting 
in a larger contrast defect at the main beam splitter. 



III. ANALYTICAL DESCRIPTION OF MODE 

COUPLING VIA MIRROR SURFACE 

DISTORTIONS 

Using Zernike polynomials as a description of mirror 
surface distortions we can look at the coupling between 
Laguerre-Gauss modes analytically. The coupling be- 
tween different Laguerre-Gauss modes when the surface 
is described by a particular Zernike polynomial is given 
by: 



p,l,p' ,1' 
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(7) 



In order to simplify the integral we use the fact that when 
kZ is small we can approximate: 

exp {2ikZ) ?^ 1 + 2ikZ 

The amplitudes of the Zernike polynomials in the mirrors 
used in gravitational wave detectors are not expected to 
exceed 10 nm. With a wavelength of 1064 nm we have 
2k Z ^0.1 and so the approximation should be suitable 
for this investigation. We are concerned with coupling 
into other modes, not back into the input mode, so the 
equation to solve becomes: 
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due to the orthogonal properties of LG modes. 

Both the Zernike polynomials and the Laguerre-Gauss 
modes can be easily separated into their angular and ra- 
dial parts. The angular function to integrate is: 
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Figure 4: Plots of the intensity distributions of the order 9 hehcal Laguerre-Gauss modes. From left to right: LGo±9, 
LG2±5, LG3±3 and LG4±i. 
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Considering the even Zernike polynomial we obtain: 
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As the integral is evaluated over the entire unit disc and 
gzo _ ^iNx27T _ -^^ where N is an integer, the integral is 
equal to 0. The only combination of Zernike polynomials 
and Laguerre-Gauss beams to give a non-zero result oc- 
curs when one of the exponentials disappears before the 
integration takes place. This occurs when we have: 
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(9) 



The same condition also gives the only non-zero results 
for the odd Zernike polynomials. This is a very interest- 
ing result as it suggests that surfaces described by Zernike 
polynomials will only cause significant coupling from one 
LG mode to another if the Zernike azimuthal index is 
equal to the difference between the azimuthal indices of 
the two modes. This requirement for m also gives the 
minimum order (n) of Zernike polynomial required to 
produce significant coupling, as r?2 < n. 

Using this condition we can integrate with respect to 
(p. The integrals were found to be tt for the even Zernike 
polynomials and ±i7r for the odd polynomials, depending 
on the sign of (/ — V). The final equation is given by: 
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where X = ^^, R is the Zernike radius and w is the 
beam radius, and 7 is the lower incomplete gamma func- 
tion. The full derivation is given in Appendix [Aj 

In our approximation of the coupling coefficients we 
only consider the magnitude of the coefficients. How- 
ever, the real coefficients have both real and imaginary 
parts indicating that there is some phase shift caused by 
the distortions. Therefore, when considering the coupling 
from a surface in terms of the coupling from the individ- 
ual polynomials making up the surface we also need to 
consider the phase shifts. The largest possible coupling 
from a surface occurs when all the individual Zernike 
couplings have the same phase and therefore the magni- 
tude of the coupling is equal to the sum of the individual 
couplings. 



IV. ANALYSIS OF COUPLING INTO ORDER 9 
MODES 



We want to verify the results of the analytical descrip- 
tion of the coupling coefficients. Using the condition for 
significant coupling outlined previously in Sec. [ill] we can 



(10) 



identify the azimuthal Zernike indices which will cause a 
large amount of coupling from LG33 into the other or- 
der 9 modes. These are summarised in table [TTl Because 
m < n this condition for the azimuthal index also tells 
us the lowest order (n) Zernike polynomial required to 
cause a large amount of coupling from the LG33 beam 
into each of the other order 9 modes. 

Higher order Zernike polynomials represent higher or- 
der spatial frequencies, which generally have smaller am- 
plitudes in the mirror surfaces. Therefore we would ex- 
pect the coupling caused by higher order polynomials, 
such as into LGi_7 and LG0-9, to be smaller than those 
caused by lower order polynomials. We would also ex- 



m 2244668 10 12 

Up^i mode 2, 5 4, 1 1, 7 4,-1 0, 9 3,-3 2,-5 1,-7 0,-9 



Table IL The azimuthal index (m) of the Zernike polynomial 
required to achieve significant coupling from an LG33 incident 
beam into the other order 9 modes. 




Figure 5: Plots of the coupling coefficients for different order 
9 LG modes when an LG33 beam is incident on a surface 
described by the Zernike polynomial Z|. The amplitude of 
the coefficients is plotted against the ratio of the beam radius, 
w, and the Zernike radius, R. These plots show numerical 
results for the coefficients without any approximation. The 
coefficients for LG4-1 and LG17 are significantly larger than 
those of the other modes. 



pect the polynomials with m = 2, 4 and 6 to have a large 
effect on the beam purity as they each couple from LG33 
into two other order 9 modes. 

Using Matlab the original integration (Eq. u\ was per- 
formed numerically, computing the coupling occurring 
from a mirror surface defined completely by a single 
Zernike polynomial. This particular example shows the 
results for Z4, in which we expect a large coupling from 
LG33 into LG17 and LG4_i (table |ll|) and much less cou- 
pling into other order 9 modes. The coupling coefficients 
between the LG33 beam and all the other order 9 modes 
were calculated. The numerical integration was carried 
out for a range of ^ and the results, for a polynomial 
amplitude of 1 nm, are summarised in Fig. [5] In this 
plot the two largest coupling coefficients over this range 
of ^ correspond to LG4_i and LG17. The coefficients 
for the other LG modes are significantly smaller. These 
results agree with our predictions. 

Fig. 6]shows a comparison of the analytical results from 
Eq. [To with the numerical results. Here only the coupling 
coefficients for LG4_i and LG17 are plotted as the analyt- 
ical approach gives for 772 7^ \^ ~^'\' Over this range the 
two sets of numbers match up very well and we consider 
this a good confirmation of the analytical approximation. 



A. Beam size and Zernike order 

Figs. [5] and |6] seem to suggest that the coupling co- 
efficients also depend on the beam size relative to the 
radius of the Zernike polynomial, or mirror radius. How- 
ever, this ratio is typically not a free parameter. For a 
good reduction in thermal noise a large beam radius is 




Figure 6: A plot showing the coupling into the LG4-1 and 
LG17 modes from an LG33 beam incident on a surface de- 
scribed by the Zernike polynomial Z|. The results are plotted 
against the relative beam radius on the mirror, ^ . The results 
from an analytical approximation and a numerical integration 
method are plotted. 



desirable. An upper limit for the beam width can be 
derived from optical loss due to beam clipping. This 
so-called clipping loss refers to the power lost over the 
mirror edges given by [7 : 



clip 
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(11) 



where the integral represents the normalised power re- 
flected by a perfect mirror of finite size (see Appendix [B|). 
In gravitational wave detectors the clipping loss should 
be lower than 100 ppm (parts per million) and often an 
arbitrary requirement of 1 ppm is used during the inter- 
ferometer design phase. This yields an optimal beam ra- 
dius of 0.232 R {R being the mirror radius) for an LG33 
beam. A clipping loss of 100 ppm instead leads to an 
optimal beam radius of 0.255 R. 



V. CAVITY SIMULATIONS WITH ADVANCED 
LIGO MIRROR MAPS 

We want to investigate how the degeneracy of the or- 
der 9 modes affects the purity of an LG33 mode in high 
finesse cavities. The aim of this investigation is to assess 
the effects of higher-order mode degeneracy and derive 
requirements for the mirror surfaces which would result 
in an acceptably high LG33 beam purity. The Advanced 
LIGO cavities consist of two curved mirrors; the input 
test mass (ITM) and the end test mass (ETM) separated 
by an arm length of approximately 4 km (Fig. I7|. The 
design properties of these mirrors are summarised in ta- 
ble III |T8l [19], giving a high finesse of 450. 

To simulate mirror surface distortions we used mir- 
ror maps measured from uncoated mirror substrates pro- 
duced for Advanced LIGO [25 . However, the Advanced 
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Figure 7: The optical layout of an Advanced LIGO arm cav- 
ity. Light from the beam-splitter is incident on the flat anti- 
reflective surface of the ITM. The light circulates between the 
two highly reflective, curved surfaces of the ITM and ETM. 
Light is transmitted by the cavity through the ETM. 



Mirror ITM ETM 
Ra 20ppm 500 ppm 
La 1 ppm 1 ppm 
Th 0.014 5 ppm 
0.3 ppm 0.3 ppm 



Lh 



Table III: Optical properties of the mirrors designed for the 
Advanced LIGO arm cavities. The two mirrors have a design 
thickness of 200 mm and index of refraction of 1.45. R, T 
and L refer to the reflectivity, transmission and loss of power 
at the mirror, a and h refer to the ant i- reflective and highly 
reflective coated surfaces of each mirror. 



LIGO cavities were not designed to be compatible with 
the LG33 mode. The LG33 mode is more spatially ex- 
tended than the LGqo mode, and so experiences a larger 
clipping loss at the mirrors for a given beam spot size 
value. For the Advanced LIGO cavity parameters, the 
clipping loss for the LG33 mode is much larger than ac- 
ceptable, at around 35%. It was therefore necessary to 
adjust the cavity length to bring the LG33 clipping to a 
level that allowed us to carry out a meaningful investiga- 
tion. Using a cavity length of 2802.9 m we achieve similar 
clipping losses to those experienced by LGqo in the orig- 
inal cavities, for the 30 cm aperture represented by the 
mirror maps. The results from this optical setup should 
be representative of longer cavities with larger mirrors. 
The simulated cavity parameters are summarised in ta- 
ble llVl 



A. Laguerre- Gauss mode purity with Advanced 
LIGO mirror maps 

For the purposes of this investigation the mirror map 
corresponding to the Advanced LIGO end test mass 
ETM08 was used; see Fig. |2] To predict the direct cou- 
plings from this mirror we look at the Zernike polynomi- 
als representing the mirror surface. This is achieved by 
performing a convolution between the surface defined by 
the mirror map and the different Zernike polynomials: 
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(12) 



where Zmap is the surface defined by the mirror map and 
A^ is the amplitude of the corresponding Zernike poly- 
nomial in the surface. The convolution was performed 
for all Zernike polynomials with n < 30. The polynomi- 
als which cause significant coupling into the other order 
9 modes (m = 2, 4,. ..12) are summarised in table |V| 
Here the polynomials are ranked in order of the power 
they couple into the other order 9 modes when an LG33 
beam is reflected from a surface described by the polyno- 
mial. From this we can suggest which order 9 LG modes 
will have significant amplitudes in the simulated cavity. 
The two polynomials which cause the largest individual 
power couplings have m = 2. The astigmatism (Z2) in 
particular extracts a large amount of power from LG33. 
Therefore, we would expect the LG41 and LG25 modes 
to have relatively large amplitudes in the cavity as the 
m = 2 polynomials cause significant coupling into these 
modes. 

The cavity defined in Sec. [V| was simulated using the 
interferometer simulation tool Finesse [26] [211 [22]. An 
input beam of pure LG33 was used with the ETM08 mir- 
ror map applied to the end mirror and a perfect input 
mirror. The cavity was tuned to be on resonance for the 
LG33 mode and the beam circulating in the cavity was 
detected. A plot of the circulating field is shown in Fig.[8| 

The purity of an LG mode, Up^i^ in a given beam U is 
defined as \cp^i\'^ where [23]: 
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(13) 



The plot in Fig. [s] (compared with the plot of LG33 in 
Fig. PI suggests that the circulating beam now contains 



Parameter ITM Re ETM Re Cavity length 
Value [m] -1934 2245 2802.9 



Table IV: Cavity parameters for simulations of Advanced 
LIGO style arm cavities [20]. The length of the cavity was re- 
duced from the original length of 3994.5 m to prevent a large 
clipping loss when using LG33 beams. 



z;r 


polynomial 


2, 2 


4, 2 


4, 4 


6, 2 


10, 8 


other 


Am 


[nm] 


0.908 0.202 


0.213 


0.124 


0.116 


- 


Power [ppm] 


4.66 


0.331 


0.0431 


0.0099 0.0059 


< 0.005 



Table V: Zernike polynomials present in the Advanced LIGO 
mirror map ETM08 which cause significant coupling from 
LG33 into the other order 9 LG modes (m = 2, 4,. . . 12). The 
power coupled from LG33 into the other order 9 modes by re- 
flection from surfaces described by the individual polynomials 
is included, calculated from a coupling approximation. 
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Figure 8: A plot of the field circulating in a simulated high 
finesse cavity. In the simulation the Advanced LIGO mirror 
map ETM08 is applied to the end test mass and a pure LG33 
beam is injected into the cavity. 



modes other than LG33. The purity of the LG33 beam in 
this simulated cavity was found to be 88.6%. The power 
in the different modes present in the circulating field are 



summarised in table VI The results of the decomposition 
show that the surface distortions of the end mirror cause 
significant coupling into other LG modes, particularly the 
other order 9 modes. The distortions also cause coupling 
into modes of other orders, but these are not resonant 
at the same cavity tuning as the LG33 mode and so are 
strongly suppressed. 

Other than LG33 the two largest modes in the cavity 
are LG41 (5.7 %) and LG25 (5.0 %). This agrees with 
the predictions made by studying the Zernike content of 
the ETM08 mirror map (table [V|). LG4-1 and LG17 also 
have relatively large amplitudes in the cavity. This may 
be due to the direct coupling out of the LG33 caused by 
the Z4 polynomial. However, the LG4_i and LG17 modes 
can also be strongly coupled out of the LG41 and LG25 
modes via m = 2 polynomials, further contributing to the 
effect these polynomials have on the purity of the beam. 
Overall the coupling process in a cavity is complicated 
by these multiple cross-couplings, but the results of this 
simulation suggest that the direct coupling from a mirror 
surface is the dominant effect on the mode content of 
the circulating beam. A theoretical understanding of the 



Up^i mode 


3, 3 


4, 1 


2, 5 


4,-1 


1, 7 


other 


msc 




- 


2 


2 


4 


4 


- 


Power 


(%) 


88.6 


5.70 


5.02 


0.333 0.313 


< 0.05 



Table VI: The power in the LG modes circulating in a cav- 
ity simulated with Finesse with an LG33 input beam and 
the Advanced LIGO mirror map ETM08. msc refers to the 
azimuthal index of the Zernike polynomial required to cause 
significant coupling from LG33 to the given mode. 



direct coupling has therefore allowed us to make valid 
predictions about the resulting mode content. 



B. Frequency splitting 

The presence of mirror surface distortions not only 
causes coupling between LG modes but introduces ad- 
ditional phase shifts of the modes. This results in slight 
shifts of the resonance frequency of individual modes. 
These shifts in resonance frequency depend on the par- 
ticular mode and so modes of the same order will be 
resonant at slightly different frequencies. Thus the mode 
degeneracy can be broken. We will refer to this effect 
as frequency splitting. For the frequency splitting to be 
effective the shifts in resonance frequency must be larger 
than the cavity bandwidth in order to separate the reso- 
nance peaks of the different order 9 modes. 

Using the ETM08 mirror map the high finesse cavity 
defined in Sec. [V| was simulated. The beam circulating in 
the cavity was detected as the laser frequency was tuned 
around the cavity resonance. The maximum power of 
the order 9 modes in the cavity and the difference in 



their resonance frequencies is summarised in table VII 



The frequency splitting is of the order of 10 Hz, smaller 
than the cavity bandwidth of 120 Hz and so is not suf- 
ficient to completely break the degeneracy. The result 
is 10 quasi- degenerate modes. The resonance frequencies 
are slightly different and so when the cavity is tuned to 
the resonance of the LG33 mode the other modes will be 
slightly suppressed. However, the frequency splitting is 
small and therefore these modes will still have relatively 
large amplitudes in the cavity. The coupling into order 9 
modes is therefore still the dominant effect on the mode 
purity. 



C. Laguerre-Gauss mode purity Avith improved 
mirror maps 

We want to find ways in which the mirror maps could 
be improved for the specific application of the LG33 
beam. By considering the mirror surfaces analytically 
it appears that reducing the astigmatism for this par- 
ticular map could improve the purity of LG33 in our 
simulated cavity. To demonstrate this the astigmatism 
was removed from the ETM08 mirror map. The cavity 
was then simulated with this processed map, with the 
resulting circulating field detected and decomposed in to 
LG modes. The circulating beam is plotted in Fig. [9] 
Simply comparing this plot with the original circulating 
beam (Fig. Isl) suggests that the purity of the field has 
increased. The LG content of the beam is summarised 
in table |VIIll The purity of the circulating LG33 beam is 



now 99.5%, a significant improvement from the original 
results. 

The results of the decomposition show that the power 
in both the LG41 and LG25 modes has decreased signif- 



Up J mode 


3, 3 


4,1 


2,5 


4, -1 


1,7 


3, -3 


0, -9 


2, -5 


0, 9 


1,-7 


Power [W] 


221.8 


14.38 


12.62 


0.865 


0.802 


0.102 


0.038 0.038 0.032 


0.001 


Frequency [Hz] 





1.3 


0.7 


7.0 


6.6 


7.6 


-7.3 


-16.1 


-14.8 


9.3 



Table VII: A summary of the power in each order 9 mode in a simulated cavity where the mirror map ETM08 is applied to the 
end mirror and the input mode is LG33. The frequency is tuned around the cavity resonance and the difference in resonance 
frequency, compared to LG33, is included for each mode. 
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Figure 9: A plot of the field circulating in a cavity simulated 
with an input mode of pure LG33 and the Advanced LIGO 
ETM08 mirror map, with astigmatism removed. 



icantly, as predicted. The power in the other modes has 
also noticeably decreased. This result suggests that the 
astigmatism is a major factor in coupling from the LG33 
mode, not only for its direct coupling into LG25 and LG41 
but for the coupling from these new modes into other 
modes of order 9. For this setup we can conclude that 
the astigmatism should be limited in the mirror surfaces 
to reduce the problems caused by higher order mode de- 
generacy. 



Up,i mode 


3, 3 


4, 1 


2, 5 


1, 7 


4,-1 


0,-9 


other 


msc 




- 


2 


2 


4 


4 


12 


- 


Power 


[%] 


99.5 0.231 


0.208 0.0524 0.0165 0.0137 


< 0.01 



Table VIII: The power in the LG modes circulating in a sim- 
ulated cavity. The cavity was simulated with an LG33 input 
beam and with the ETM08 mirror map with astigmatism re- 
moved, msc refers to the azimuthal index of the Zernike 
polynomial required to cause significant coupling from LG33 
to the given mode. 



D. Mirror requirements for LG33 

In order to use LG33 in GW detectors we require cer- 
tain Zernike polynomials in the mirrors to be smaller 
than in the current state of the art mirrors, in order to 
achieve an acceptable beam purity in the cavity. Here we 
investigate the direct coupling from ETM08 and suggest 
limits to the amplitudes of specific polynomials in the 
mirror surfaces, presenting an Advanced LIGO mirror 
map adapted for the use of LG33. 

Using our theoretical analysis we can identify the par- 
ticular Zernike polynomials to reduce in the mirrors. We 
have already seen that the Zernike polynomials with odd 
values of n and with m > 12 don't have a large effect 
on the purity. To assess the other polynomials we use 
Eq. [To] to approximate the coupling into order 9 modes 
caused by the polynomials in the ETM08 mirror map, for 
the optical setup defined in Sec. [V| For each order 9 LG 
mode the coupling was calculated for the Zernike poly- 
nomials with n = 2, 4 ... 30 and with m required to give 
significant coupling. Fig. [lO] represents these coupling 
coefficients. This chart shows that the largest couplings 
occur from Z2, into LG25 and LG41, as previously sug- 
gested. There is also some strong coupling from Z4 and 
Z4. The other couplings are significantly smaller. There- 
fore, the first step in modifying the mirrors for LG33 is 
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Figure 10: A bar chart showing the coupling, k^'^ , ^/ , into the 
order 9 modes when an LG33 beam is incident on a surface 
described by the Zernike polynomial with m required to cause 
significant coupling. The Zernike amplitudes correspond to 
those in the ETM08 mirror map. 




Figure 11: Plots showing an approximation of the couphng 
from LG33 into the other order 9 LG modes from the ETM08 
mirror map as higher order (n) Zernike polynomials are in- 
cluded in our model. The sum of the coupling from polyno- 
mials with order < n are plotted for each mode. 



to limit these 3 polynomials to give similar couplings to 
the higher order polynomials. 

Another plot illustrating the coupling from this map is 
shown in Fig. [TT] In this plot the maximum possible cou- 
pling from the ETM08 mirror map is estimated using our 
analytical approximation. For each order 9 mode the sum 
of the coupling is calculated for all Zernike polynomial or- 
ders smaller than n. The plot illustrates the magnitude 
of the direct coupling expected as we include higher order 



Zernike polynomials in our model. Fig. 11 can be used 
to illustrate how this particular map, ETM08, should be 
adapted for LG33. The coupling into LG25 and LG41 is 
around 10 times greater than any other coupling. The 
coupling into these two modes is also not significantly in- 
creased by including the higher order modes. Therefore, 
limiting the lower order polynomials with m = 2 will 
greatly reduce the coupling into these two modes and 
the overall coupling into order 9 modes. From this plot 
we conclude that the overall coupling into order 9 modes 
can be reduced by around a factor of 10 by reducing the 
lower order polynomials. Reducing the coupling further 
will involve limits on multiple polynomials. 

By assessing the direct coupling from ETM08 we can 
set requirments for the lower order Zernike polynomials. 
When an LG33 beam is incident on the ETM08 mirror 
map in our setup the reflected beam is predominantly 
LG33. However, we find 31 ppm (parts per million) of 
the power is now in other modes, with 6.8 ppm in the 
other order 9 modes. Table |V| shows the power cou- 
pled from surfaces described by the polynomials present 
in the ETM08 map, from LG33 into the other order 9 
modes. For this map we consider the polynomials caus- 
ing a large amount of coupling as those coupling more 
than 0.01 ppm into the order 9 modes; Z2, Z4 and Z4. For 
LG33 we require these polynomials to be limited in the 
mirror surfaces. The requirements for these polynomials 
were calculated to give power couplings of 0.01 ppm into 



Z;r polynomial 2, 2 4, 2 4, 4 
Amplitude [nm] 0.042 0.035 0.100 



Table IX: A summary of the amplitude requirements for the 
Zernike polynomials required to give individual couplings of 
0.01 ppm from LG33 into the other order 9 modes in the 
ETM08 mirror map. 



the other order 9 modes and are summarised in table lIXl 
These amplitude limits were applied to the ETM08 mir- 
ror map, resulting in coupling of 19 ppm into modes other 
than LG33 and 0.043 ppm into the other order 9 modes. 
The cavity defined in Sec. [V| was simulated with this lim- 
ited map, resulting in 815 ppm impurity in the circulating 
beam. This is a very good improvement on the original 
impurity of 0.114, illustrating that a high beam purity is 
achievable with these mirror requirements. To achieve an 
even higher beam purity will involve reducing the ampli- 
tudes of these polynomials further, as well as additional 
Zernike requirements. 



VI. CONCLUSION 

We have investigated the coupling which occurs when 
Laguerre- Gauss modes are incident on a mirror with sur- 
face distortions. Taking an analytical approach we used 
Zernike polynomials to represent mirror surface distor- 
tions and derived an approximate equation for the sig- 
nificant coupling when an LG mode is reflected from a 
surface defined by a particular Zernike polynomial. This 
derivation resulted in a condition for significant coupling, 
m = \l — l'\^ where m is the azimuthal index of the Zernike 
polynomial and / and I' are the azimuthal indices of the 
incident and coupled modes respectively. This is a signif- 
icant result as it allows us to predict which order 9 modes 
will be largely coupled by particular Zernike polynomials 
and suggest which modes will have large amplitudes in 
the arm cavities 

We investigated the performance of LG33 in high fi- 
nesse cavities by simulation with Advanced LIGO mirror 
maps. This illustrated the degraded purity of the circu- 
lating beam in realistic cavities due to higher order mode 
degeneracy. The results were then analysed by looking at 
the Zernike polynomials representing our example mirror 
map. The analysis and results were consistent with the 



predictions made from Eq. [10] This suggested that astig- 
matism was causing a significant amount of coupling, 
particularly into the LG41 and LG25 niodes. This was 
confirmed when the cavity was simulated again with the 
astigmatism removed from the mirror map and we ob- 
served a dramatic increase in the LG33 mode purity. 

The analytical description enabled us to identify the 
specific Zernike polynomials which cause large couplings 
as well as the LG modes which would dominate as a re- 
sult. Using this we were able to derive certain require- 
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merits for our example mirror map, ETM08, in terms 
of limits on the amp litud es of the Zernike polynomials 
Z2, Z4 and Z4 (table IX). Using this map the resulting 
circulating beam impurity was found to be 815 ppm, a 
significant reduction from the original impurity of 0.114. 



This investigation has demonstrated that a high beam 
purity is achievable using an LG33 beam when modifica- 
tions are made to the low order Zernike polynomials in 
Advanced LIGO mirrors. Implementing the LG33 beam 
in gravitational wave detectors will be challenging as we 
require very small amplitudes on these lower order poly- 
nomials. We should also consider that the example mir- 
ror surfaces considered here refer to uncoated substrates. 
The coating process is likely to add to the lower order 
features in the mirror surfaces. However, using this an- 
alytical approach we can derive specific requirements for 
the mirror surfaces leading to designs for suitable mirrors 
for these higher order beams. 
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Appendix A: Derivation of coupling coefficients 



The product of two Laguerre- Gauss modes is: 



Up,iUpf^if = 



1 2 

w^7r\lj\l\+p)KW\+pOl 
xexp(i(2p+|/|-2p^-|r|)^) 



\/2r 
w 



\i\ + \i'\ 






(Al) 



exp 



2r' 



exp (i(/) (/ — l')) 



The following derivation follows from Eq. [9J Currently 
we are concerned with the magnitude of the coupling 
coefficients, so we ignore any constant phase shifts and 



integrate with respect to (j): 
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(A2) 



In order to further simplify the equation the following 
variable substitution is made: 



2r' 

and a new limit to the integral: 

2R^ 



X 



w^ 



(A3) 



(A4) 



where R is the Zernike radius. This gives the integral: 







pip' I 



+p)\(\i'\+p')\ 
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(A5) 



X L^p^{x)L^p,^{x)exp{-x)dx 



Substituting in the sums representing the Laguerre poly- 
nomials and the radial Zernike function as a function of 
X gives: 

KZ',i' = A^kVp\p'\(p+\i\)Kp' + \i'\y. 



P p' \{n-m) 

1 

{p^-j)\m^j)\j\a{n^m)-h)\ (^^) 



{p - i)\{\l\ ^ i)\i\ 
{n-h)\ 



(^(n-m)-/i)!/i!xi(^-2^) 

/''x^+^-^+^(l^l+l^'l+^)exp(-x)dx 
Jo 



This type of integral results in the lower incomplete 
e;amma function: 



7(a,x) = / e-^e-^dt 
Jo 



(A7) 



When a is equal to n, an integer, the function is given by 
the following sum: 



7(n,x) = (n-l)! l-e--^ 



n— 1 iu ^ 



fe=0 



k\ 



(A8) 
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Therefore, the final equation for this approximation of 
the magnitude of the coupling coefficients is given by: 



p p' \{n-m) 
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EE E 

1 
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(A9) 



Appendix B: Clipping loss 

Clipping loss is given by: 



1 



clip = -'^ / I'^P'^I 



|t/„ 



(Bl) 



where the integral represents the normalised power re- 
flected by a perfect mirror. S defines an infinite plane 
perpendicular to the beam axis. The magnitude squared 
of an LG mode is: 



\U^ 



pM 



1 2p\ f2/' 

w'^ 7r(|/| -\- p)\ \ w 



2\ \i\ 



"•^''(S)) ^^P(-S 



(B2) 



where w is the beam radius, p and / are the mode in- 
dices and r is the radial position. Integrating over the 
surface and taking a similar approach as for the coupling 
coefficients we get: 



«ciip = i-pKp+i^i)!EE 



(-1) 



n-\-'m 



m—0 n=0 



{p — n)\{p — m)\ 



+ n)!(|/| ^m)\n\m\ 



7(|/|+n + m + l,X) 
(B3) 



where X = ^^, R is the radius of the mirror and 7 is 
the lower incomplete gamma function. 



Appendix C: Zernike composition of a mirror surface 

For certain Zernike polynomials (those with non-zero 
m) their amplitudes in a surface depend on the orienta- 
tion of that surface with respect to the Zernike surface. 
For example, consider the two polynomials responsible 
for astigmatism, Z2 • The two polynomials actually de- 
scribe the same shape, with one just rotated by 90 degrees 



with respect to the other. Therefore, rotating a surface, 
such as the one described by mirror map ETM08, will 
change the amplitudes of these two polynomials within 
the surface. Fig. [12] illustrates this effect. The plots 
show the amplitudes of the order 2 Zernike polynomi- 
als present in the ETM08 mirror surface as it is rotated. 
As expected the amplitude of the Z2 polynomial remains 
constant as it has no angular dependence. The ampli- 
tudes of the Zf^ polynomials oscillate and, at a certain 
rotation (around 120°) the astigmatism of the surface is 
completely described by Z2"^, and 90° later completely 
described by Z^ . The root mean squared amplitude 
(A2) of the polynomials remains constant. 
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Figure 12: Plots of the amplitudes of the order 2 Zernike 
polynomials present in the surface described by the ETM08 
Advanced LIGO mirror map as a function of the rotation of 
the surface. 



Appendix D: Analysis of Zernike approach 

We have used Zernike polynomials to describe mirror 
surface distortions and analyse the coupling that occurs 
from LG33 into other order 9 modes. This approach ap- 
pears to be very suitable as we have been able to identify 
specific polynomials which extract significant amounts of 
power from the input mode. However, there is an alter- 
native method used to investigate mirror surface distor- 
tions, which involves looking at the spatial frequencies 
present in real mirrors. In this section we compare these 
two methods. 

To look at the spatial frequencies present in realistic 
mirrors we perform a 2D Fourier transform of the sur- 
face height data. The resulting spectra is then analysed 
and synthetic maps are created with the same spatial fre- 
quencies. This method focuses on identifying particular 
spatial wavelengths which cause a large degree of cou- 
pling from LG33. Many synthetic maps are created and 
used in simulations of gravitational wave detectors. A 
statistical approach is then taken to determine the ex- 
tent of the coupling when specific spatial frequencies are 
present in the mirror. 
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In the Zernike approach we look at the different poly- 
nomials present in mirror surface distortions. This can be 
thought of as equivalent to looking at the spectra of the 
mirror surfaces as the different polynomials represent dif- 



interesting polynomials and modes. 



ferent spatial frequencies. The plot in Fig. [13] illustrates 
this. The spectrum of the ETM08 Advanced LIGO mir- 
ror map is shown, along with the spectra of maps made 
up from the Zernike polynomials present in the ETM08 
mirror. Each of the Zernike maps recreates the LIGO 
map with polynomials up to a certain order. The plots 
show that as the order of Zernike polynomials present in- 
creases the higher spatial frequencies are represented in 
the mirror map. This is because these higher order poly- 
nomials represent the higher order spatial frequencies. 
Looking at the spatial frequencies present in the Zernike 
polynomials we found that the frequencies depended on 
the order, n. A consequence of this is that if we just con- 
sider the spatial frequencies present in the mirror maps 
we will not be able to distinguish between polynomials 
with different m. As we have seen, the azimuthal in- 
dex is very significant as it determines which modes are 
largely coupled from LG33. Therefore looking at the spa- 
tial frequencies doesn't identify the important shapes in 
the mirror surfaces. The Zernike approach would seem 
to be the most suitable as this allows us to identify the 
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Figure 13: A plot showing the spectrum of the Advanced 
LIGO mirror map, ETM08 and the spectra of maps created 
from Zernike polynomials present in ETM08. The Zernike 
maps go up to a certain maximum order, recreating higher 
spatial frequencies with higher orders. 
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At the time of our analysis the coated mirrors were not 
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